Probing the Primordial Power Spectrum with Cluster Number Counts 
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We investigate how well galaxy cluster number counts can constrain the primordial power spec- 
trum. Measurements of the primary anisotropies in the cosmic microwave background (CMB) may 
be limited, by the presence of foregrounds from secondary sources, to probing the primordial power 
spectrum at wave numbers less than about 0.30 hMpc" 1 . We break up the primordial power spec- 
trum into a number of nodes and interpolate linearly between each node. This allows us to show 
that cluster number counts could then extend the constraints on the form of the primordial power 
spectrum up to wave numbers of about 0.45 ftMpc . We estimate combinations of constraints from 
PLANCK and SPT primary CMB and their respective SZ surveys. We find that their constraining 
ability is limited by uncertainties in the mass scaling relations. We also estimate the constraint from 
clusters detected from a SNAP like gravitational lensing survey. As there is an unambiguous and 
simple relationship between the filtered shear of the lensing survey and the cluster mass, it may be 
possible to obtain much tighter constraints on the primordial power spectrum in this case. 



I. INTRODUCTION 



A crucial element of cosmology is the form of the 
primordial fluctuations. These fluctuations provide the 
seeds for structure formation which we observe to- 
day through the cosmic microwave background (CMB), 
galaxy, and cluster surveys. Current data is consistent 
with the primordial fluctuations being scalar, adiabatic, 
Gaussian and having a power spectrum with a simple 
power law parameterization [l], [2j • The primordial fluctu- 
ations may have been generated during a period known as 
'inflation', where the accelerated expansion of the primor- 
dial Universe is driven by a potential dominated scalar 
field or fields (see for example Liddle and Lyth Q). If 
inflation was driven by a single scalar field with a smooth 
potential, then the power spectrum of primordial fluctu- 
ations is predicted to be generally quite close to a power 
law form, although in some cases there may be significant 
running of the spectral index. However, if inflation was 
driven by multiple fields or by a single field with a fea- 
ture in its potential, then the primordial power spectrum 
may contain hills, valleys, oscillations or other features 
(see for example 0, H, IE 0, H H ) • 

The two main approaches to probing the primordial 
power spectrum are either to assume a specific form 
for a feature in the primordial power spectrum (see for 
example [E EE El, 113, 113) or to try and reconstruct 
the primordial power spectrum non-parametrically (see 
for example 0, EI EE EE EE EE El)- The cleanest 
probe of the primordial power-spectrum is the cosmic mi- 
crowave background (CMB). However, it is probably lim- 
ited to wave numbers smaller than about 0.30 /iMpc -1 
as beyond that foreground contamination from secondary 
sources arc likely to dominate the cosmological signal. 

In this article, we investigate to what extent galaxy 
cluster number counts can probe the primordial power 
spectrum. Traditionally, the effect of the primordial 
power spectrum on clusters is probed with the constraints 
on erg , which is the dispersion of the linear theory matter 
fluctuations smoothed on scales of 8 /i -1 Mpc. However, 



as we will discuss in Sec. Ill Al the value of ag encompasses 
a very broad range of wave-numbers and so it is advan- 
tageous to break up the primordial power spectrum into 
effectively several small bins. Additionally, as is sensi- 
tive to other cosmological parameters such as the matter 
density, the dark energy equation of state, primordial 
non-Gaussianity, and non-zero neutrino mass. Rather 
than simply comparing a$ inferred from clusters with 
that from inferred from the CMB, it may be better to 
introduce new parameters to account for the possible de- 
viations from the fiducial model of ACDM consisting of 
a featureless, adiabatic, and Gaussian primordial power 
spectrum. In this article we investigate how well the com- 
bination of the CMB and cluster surveys can constrain 
the primordial power spectrum to be featureless. 

We begin in Section |TT] with contrasting how the clus- 
ter number counts and the primary CMB probe the pri- 
mordial power spectrum. In Section IIIII we forecast the 
constraints on the primordial power spectrum from the 
Sunyaev-Zeldovich effect (SZ) PLANCK[Ii[ and South 
Pole Telescope [42|] (SPT) cluster surveys and combine 
them with the forecasted constraints on the primordial 
power spectrum from the PLANCK and SPT primary 
CMB survey. We also estimate the constraints from the 
SNAP lensing cluster survey. Concluding remarks are 
given in Section ITVl 



II. 



DEPENDENCE ON PRIMORDIAL POWER 
SPECTRUM 



The dimcnsionless primordial power spectrum, as func- 
tion of the comoving wave number k, is usually parame- 
terized as a power law 



V n ,o(k) = A 2 



^-pivot 



(1) 



where the amplitude (Afj), spectral index (n s ), and pivot 
point (fcpivot) are taken to be independent of k. We model 
deviations from this form in a similar way to Bridle et al. 
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FIG. 1: The effect on the primordial power spectrum of in- 
dividually changing nodes 1 to 15 to a value of 1.1 in our 
feature function F. Each changed node is a plotted as a sep- 
arate colour and odd numbered nodes are plotted as dashed 
lines. 
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FIG. 2: Effect of individually setting the node at F[kf] = 1.1 
on the variance of the matter field smoothed with 8/i -1 Mpc 
top hat window function. 



[3 and Spergel et al. Q: 
Vn(k) = F(k)VnAk) 



(2) 



where the "feature function", F(k), is specified by the 
values at 17 nodes logarithmically spaced in fc-space with 
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(3) 



Linear interpolation in log(fc) is used to determine F be- 
tween each node. The effect of individually changing 
nodes 1 to 15 is shown in Fig. [TJ Nodes and 16 will 
always be fixed to 1. 

Spergel et al. [J (see their Fig. 11) found that WMAP 
III data constrained the nodes best at around k ~ 
0.01 h Mpc -1 with a 1 sigma error on F of about 0.3. 
While at k ~ 0.1 /iMpc -1 the WMAP III data provided 
practically no constraint due to its relatively large beam 
size and measurement noise. It would be interesting to 
see how these constraints could be improved with the ad- 
dition of small scale CMB experiments and other large 
scale structure data such as the Sloan Digital Sky Survey, 
but that is beyond the scope of the present paper. It is 
probably reasonable to say that F is not known to better 
than 10% precision with current data, certainly for larger 
wave numbers (k > 0.1 /iMpc -1 ) where the current con- 
straints are likely to be significantly more uncertain. For 
this reason we take F(ki) — 1.1 to conservatively illus- 
trate the effect of node changes in Figures [I] [2l [3l and 

SI 



A. Effect on number counts 

The linear theory matter power spectrum at a rcdshift 
z is given by 



P(k,z) (xT 2 (k,z)kViz{k) 



(4) 



where T(k, z) is the matter transfer function. We use 
CAMBf2a||44t to evaluate P(k,z) and we modified the 



"ScalarPower" function in CAMB to include our feature 
function F{k). The accuracy and sample boost parame- 
ters in the CAMB initialization file where all set to the 
value of 3 which increases the precision and reduces the 
amount of interpolation. Although, there is some in- 
terpolation in k used by CAMB, the sampling is much 
smaller than the width between our nodes. We include 
lensing of the CMB, but its effect is negligible for our 
primary region of interest, I < 2000. Throughout this 
paper we assume a flat ACDM Universe with no tensor 
perturbations and we use the WMAP5 maximum likeli- 
hood parameters [2l[ : 

fl b h 2 = 0.0227, n c h 2 = 0.108, n = 0.961, r = 0.089, 

A|, = 2.41 x 10 -9 , h = 0.724, (5) 

where we have used the WMAP chosen pivot of fc p i vot = 
(l/500)Mpc -1 and the parameters have their usual 
meaning. 

The variance of the linear theory matter field, which 
has been smoothed by a top hat filter on a comoving 
length scale R, is given by 



a 2 (R,z) 



where 



W{kR) = 3 



dk k 3 
~~k 2^2 



P(k,z)W 2 (kR) 



s'm(kR) cos(kR) 



(fci?) 3 {kRf 



(6) 



(7) 



The top hat smoothing suppresses the contribution of any 
fluctuations located at wave number kf 1/R. Fig. [5] 
illustrates how changing the primordial power spectrum 
alters a when R = 8/i _1 Mpc. As can be seen, the main 
contribution is from the node at k = 0.20 h Mpc - , but 
that there is also a reasonable contribution from sev- 
eral of the neighboring nodes. As will be discussed in 
Sec. IIIB| the CMB is able to well constrain the primor- 
dial power spectrum at least for k < 0.30 h Mpc - . So 
it follows that as is mainly affected by scales which are 
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constrainable by the CMB. That is why our method of 
breaking up the primordial power spectrum into linearly 
interpolated nodes is useful as it gives enough flexibility 
to separate the parts of the primordial power spectrum 
which affect cluster counts but may not be constrainable 
by the CMB. Also, as is sensitive to not only the pri- 
mordial power spectrum but also the other cosmological 
parameters, such as Q m = f2 c + fib, the dark energy equa- 
tion of state, primordial non-Gaussianity, and non-zero 
neutrino mass. 

The number density of halos (bound objects) may be 
predicted using the smoothed linear theory density field 
|22| ■ For a background non-relativistic matter density of 
Pm = fimPtotai, the number density (n) of halos of mass 
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depends, to a good approximation, on the primor- 
dial power spectrum only through its effect on a(R, z) 



dn(z) Pm dlncr(z)" 
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We use the Sheth-Tormen mass function [22 
which 
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for 



(10) 



where A = 0.3222, a = 0.707, p = 0.3, and S c = 1.686. 
The top hat smoothing in Eq. ^ suppresses the con- 
tribution of any change to the primordial power spec- 
trum located at wave number kf ^> 1/R. Combined 
with Eq. I|8|). this implies that a change in the primor- 
dial power spectrum at kf has a suppressed effect on the 
number density on mass scales satisfying 



M 



h- l M, 



> 10 1 



( k 



f 







\hMpc~ 



-3 



(11) 



The number of clusters per redshift interval above some 
mass threshold M m - ln is given by 



-^W>M mia ) = f sky -U 



Lj m ^ {m -' ] - 



(12) 



Where / s k y is the fraction of the sky being observed and 
the volume element is given by 



dV 

dz 



4tt 



dz' 



n 2 



o H(z')_ 
and H{z) is the Hubble parameter 

H(z) = HoV(n m (i + z) 3 + (i-n rn )) 



(13) 



(14) 




FIG. 3: Effect of change in primordial power spectrum on the 
SNAP lensing cluster counts (M min = 1O 14 /i -1 M , / sky = 
0.024). The plots are for the featureless (F = 1) power 
spectrum (green, solid), _F[0.45] = 1.1 (blue dashed), and 
F[0.30] = 1.1 (red dotted). The vertical error bars are one 
sigma and the horizontal error bars indicate the bin width. 



The effect of a change in the primordial power spec- 
trum on the number counts is illustrated in Fig. [3] for 
a SNAP-like gravitational lensing cluster survey. The 
SNAP selection function can be approximated by setting 
M min = 10 14 h- 1 M Q in Eq. O (see Section[TTD)- Due to 
the exponential suppression in Eq. (| 10|) . dN/dz mainly 
depends on the mass scale M m ; n . From Eq. (jTTj) this 
implies the SNAP cluster survey will become insensitive 
to the primordial power spectrum for k S> 0.2/iMpc -1 . 
This is consistent with F[0.45] = 1.1 having less of an 
effect than F[0.30] = 1.1 as illustrated in Fig. [3] 



B. Effect on the CMB 

The primordial power spectrum is probed over a wide 
range of wave numbers by measurements of the primary 
CMB anisotropies (see for example Hu and Okamoto 
[l5|h Both the temperature (T) and E-mode of the po- 
larization (E) probe scalar perturbations. 



2tt 



xx' 



= I dhx kT^(k)Tf(k)V n (k) 



(15) 



where X, X 1 € T, E. The projection of a mode of wave- 
number k on to the surface of last scattering (a sphere 
of comoving radius £>») results in the CMB transfer 
functions having the form Tf- ~ je(kD^). Where ji is 
the spherical Bessel function of order £ which peaks at 
£ w kD*. Therefore a feature in the primordial power 
spectrum at wave- number kf is mapped onto a feature 
in CMB angular power spectrum at 



fc/Z?* 



10 4 



hMpc 



(16) 



We also used our modified version of CAMB to evalu- 
ate Cf x . Although, there is some interpolation in £ 
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FIG. 4: Effect of a change in primordial power spectrum on 
the CMB (TT) angular power spectrum. Error bars (yellow) 
include instrument noise, beam size, and cosmic variance ef- 
fects. They are centered on the featureless (F — 1) (black) 
power spectrum and are for SPT. The red dashed curves are 
for F[0.21] = 1.1 (top panel), F[0.30] = 1.1 (middle panel), 
and .F[0.45] = 1.1 (bottom panel). Perfect subtraction of sec- 
ondary foreground sources is assumed for the figure 



used by CAMB, the sampling is much smaller than the 
width between our nodes. Fig. Q] shows the effect F ^ 1 
on Cj T with predicted SPT error bars (see Section . 
The panels in the figure are consistent with Eqs. and 
(p~6|) in that having F{ki) = 1.1 translates roughly into 
a triangular perturbation of about 10% amplitude in Cf.. 
In Fig. [3] we have not included the foreground contribu- 
tion from secondary sources which will probably be hard 
to completely remove for I > 2000. For this reason, as 
done by Hu and Okamoto [IJ| and Leach [l||, we will 
restrict ourselves to £ < 2000 when evaluating the fore- 
casted marginalized errors. 



III. FORECASTS 

We use the Fisher matrix formalism to make forecasts 
on how well the primordial power spectrum can be con- 
strained. We take as our fiducial model the WMAP5 
maximum likelihood parameters, Eq. ([5]). We consider 
two cluster count SZ experiments: PLANCK and SPT. 
For PLANCK, the selection function can be approxi- 
mated by / sk y = 0.8 and M m in,PLANCK = 5 x lO 14 /i _1 M 
[27l ]. A common approximation for the SPT selection 
function is / s k y = 0.1 and Af m in,SPT = 1-75 x \Q 1A h~ 1 M Q 
(see for example Sefusatti et al. (28[, Lo Verde et al. 
[29(). This may be overly optimistic, but increasing it 
to Mnun,SPT = 3 x 1O 14 /i _1 M docs not qualitatively 
change our conclusions as the uncertainty in the mass 
scaling relationship turns out to be the limiting factor. 
We take z £ [0,3] with bin sizes Sz = 0.1, although the 
bulk of the constraining power comes from z < 1 and 
the results are negligibly changed if we take z < 2. We 
also forecast the SNAP cluster lensing survey constraints 
[30l ] (also see Hamana et al. [3l|, Wang et al. [32| . Fang 
and Haiman Takada and Bridle (34|). Here we take 
z < 1.5, / sky = 0.024, and M min , SNAP = W^H^Mq. 
This roughly matches the number of clusters found using 
a more accurate selection function of Marian and Bern- 
stein (30j when we use their fiducial model cosmological 
parameters. For our fiducial model, Eq. (O, we predict 
a total of 8888 clusters detected from a SNAP cluster 
survey. Also, [H[ found that the signal to noise ratio of 
a more realistic selection function was about the same as 
taking a mass limit of 1O 14 /i _1 M0. Our SNAP selection 
function is biased to slightly higher redshifts than that 
of Marian and Bernstein [301 ] , but we expect this not to 
alter our predicted constraints significantly. 

The number of clusters in redshift bin i is Poisson dis- 
tributed with the expected number (e^) given by inte- 
grating Eq. (fl2|) over the red shift bin. Element (j, k) of 
the Fisher matrix for a cluster count experiment is given 
by HI 

where pj consists of the cosmological parameters in 
Eq. ([5]) except for and n as they will be almost com- 
pletely degenerate with the feature function F. To ac- 
count for uncertainties in the mass of the SZ observed 
clusters, we also allow M m i ni pLANCK and M m i nj spT to be 
free parameters and give them both priors of 10% one 
sigma errors. The lensing observed clusters have a well 
determined mass-scaling relation and so we do not take 
.snap to be a free parameter [3l|, [H, [H, |34[. Ad- 
ditionally, we allow the nodes 1 to 15 in Eq. of F 
to vary. The derivatives are taking at the fiducial values 
of the parameters which in the case of the M m i n are the 
previously specified values and for the feature function, 
F{ki) = 1 for all i. The derivatives are approximated by 
the symmetrized form of a difference equation so as to 
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TABLE I: PLANCK Instrument Characteristics 



Center Frequency (GHz) 


70 


100 


143 


217 


6 (FWHM arcmin) 


14 


10 


7.1 


5.0 


or (a«K) 


12.8 


6.8 


6.0 


13.1 


a E (a»K) 


18.2 


10.9 


11.4 


26.7 



minimize truncation error (see for example Press et al. 
The CMB Fisher matrix is given by (see for example 

0) 



t x,x> y% FJ 



(18) 



where the covariance matrix can be obtained from Zal- 
darriaga et al. (37| and it depends on the temperature 
noise per pixel (or), the polarization noise per pixel (erg), 
the pixel area in radians squared (9 2 = 47r/A^ p i X ), and the 
beam window function which we approximate as Gaus- 
sian (Be ~ exp(— 1(1 + l)cr 2 ). The values we use are 
taken from the PLANCK blue book[45| and are listed in 
Table |T] (note that 8 needs to be converted to radians). 
We use crfc = 0/y8 log[2] and combine the different fre- 
quency bands as specified in Bond et al. [38|. We also 
include the constraints from SPT primary CMB temper- 
ature measurement for which we just use one band with 
6 = 1 arcmin and ctt = 10/iK. We take the range in I to 
be 2 to 2000. At higher £, secondary sources of temper- 
ature and polarization will likely prohibit the extraction 
of cosmological information from the primary CMB. 

The expected covariance matrix of the parameter er- 
rors is approximated by the inverse of the Fisher matrix. 
The expected marginalized one sigma error bars are then 
given by the square roots of the diagonal elements of the 
expected covariance matrix. Also, experiments can be 
combined by adding the Fisher matrices. When we com- 
bine PLANCK and SPT, primary CMB or SZ cluster 
detection, we reduce the / s k y for PLANCK to 0.7 so as 
not to count the same clusters twice. We plot the ex- 
pected one sigma marginalized errors for each node of F 
in Fig. For each node, marginalization is done over 
all other nodes, cosmological and mass parameters. The 
lack of constraint at low k is from cosmic variance. At 
high k, the primary CMB does not constrain the pri- 
mordial power spectrum as it is assumed to be limited 
to I < 2000 due to secondary sources. As shown in 
Fig. [5l including the lensing detected clusters makes a big 
improvement in constraining the primordial power spec- 
trum at k = 0.45/iMpc- 1 . The PLANCK (TT, TE, EE) 
data only constrains the primordial power spectrum at 
k = 0.45/iMpc -1 to about 250%. The constraints on the 
other cosmological parameters we included are not sig- 
nificantly altered by the addition of clusters. Combining 
all the surveys, we considered, improves the constraint 
at k = 0.45ft.Mpc _1 to about 25%. 

The main reason why the SZ cluster surveys are not as 



2.5 
2.0 
1.5 
1.0 
0.5 



10" 



10" 



10" 



10" 



10" 



10" 1 



k/h Mpc 



FIG. 5: Expected marginalized one sigma errors for the 
different nodes of feature function (F). Errors are shown 
for PLANCK (TT, TE, EE) (black triangles), PLANCK 
(TT, TE, EE) and SPT (TT) (red diamonds), PLANCK 
(TT,TE,EE, clusters) and SPT (TT, clusters) (blue squares), 
and PLANCK (TT, TE, EE) and SNAP (clusters) (green cir- 
cles). 



effective at constraining the primordial power spectrum 
as the cluster lensing survey is due to a degeneracy with 
the uncertainty in their mass parameters, see the top 
panel of Fig. [6l As the primary CMB is taken to be 
limited to I < 2000, it is barely altered by changes to the 
node at F(0Ab), (see Fig. [4j. As can be seen from the 
bottom panel of Fig. [6j the addition of the SNAP lensing 
cluster survey dramatically sharpens the constraint at 
F(0.45). This is consistent with the large effect seen in 
Fig. [3l The primary CMB would need to be removed of 
foregrounds to a high level of accuracy up to £ ~ 3000 in 
order to accurately measure F(0Ab) without the aid of 
clusters (see Fig. 0]). 



IV. CONCLUSIONS 

In this article we have investigated what role clus- 
ter number counts can play in constraining the primor- 
dial power spectrum. We found that if PLANCK and 
SPT primary CMB anisotropy measurements arc lim- 
ited to I < 2000 by secondary source foregrounds, then 
they can only measure the primordial power spectrum 
at k = 0.45 /iMpc -1 to about 220% precision due to 
the degeneracy between changes in the primordial power 
spectrum at smaller k. Including SPT and PLANCK SZ 
cluster surveys increases the precision to about 124%, but 
they are limited by a degeneracy with the determination 
of the observed clusters' masses. While a SNAP like grav- 
itational lensing cluster survey combined with PLANCK 
primary CMB data may be able to increase the precision 
to about 30% due to the accurate relationship between 
the observed lensing shear and cluster masses. 

In the current article we have used simple minimum 
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FIG. 6: Marginalized probability contours containing 95% 
of the expected posterior probability. The top panel is for 
PLANCK (TT, TE, EE) and SPT (clusters). The bottom 
panel is for PLANCK (TT, TE, EE) alone (blue solid) and 
with SNAP (clusters) (red dashed)). 

mass selection functions for the cluster surveys. As we 
are investigating a new use for cluster surveys, we think 
it is justifiable to initially get a more qualitative and 



easily reproducible forecasted constraint. In future work, 
building on the current investigation, we will evaluate 
the effect on our current conclusions of more realistic 
selection functions such as those in [3(J H^] . 

It is common in cluster surveys or forecasts of surveys 
to evaluate the constraints on as with and without the 
CMB. If there is a feature at high k in the primordial 
power spectrum then this could lead to a discrepancy in 
the value of as obtained from the CMB. Our method 
could then be used to determine the size of the feature 
needed to explain such a discrepancy. Other possible 
sources in a discrepancy between the CMB and cluster 
constraints on as could be a non-cosmological constant 
source of dark energy, primordial non-Gausianity, and 
sufficiently large neutrino mass. Degeneracies between 
the dark energy equation of state and non-Gaussianity 
where looked at by Sefusatti et al. [28j|. They found that, 
provided redshift information was available, there was 
not significant degeneracy between the two. Our method 
could be useful in determining what the observational 
degeneracies between features and other possible sources 
of discrepancy in as are. Comparing our Fig. [3] with 
Fig. 2 of Sefusatti et al. [28[ indicates that there may 
be some degeneracy between a feature in the primordial 
power spectrum and primordial non-Gaussianity. Also, 
inflation models which generate features in the primor- 
dial power spectrum may naturally generate scale depen- 
dent non-Guassianity [4(| . This could increase the overall 
change in cluster counts and thus make the feature more 
easily detectable. 

An alternative way of parameterizing features in the 
primordial power spectrum is to allow a running of the 
spectral index, dn s /dlnfc (see for example 0]). This is 
less flexible than our current approach, but may be more 
natural to implement in an inflation model. We plan to 
investigate how well cluster number counts improve the 
running of the spectral index in future work. 
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